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У статті пропонується просторова математична модель динаміки механічної системи класу одноступінчастих, 
циліндричних, косозубих, евольвентних зубчастих передач. Пропонована модель являє собою систему 
чотирнадцяти диференціальних рівнянь другого порядку відносно узагальнених координат, складену 
на основі рівнянь Лагранжа другого роду, що описують просторові коливання зубчастих коліс механічної 
системи. 


Постановка проблеми. Поглиблене дослідження динаміки зубчастих передач 
передбачає розгляд впливу різноманітних факторів на динамічну навантаженість і ко- 
ливальний процес у рамках внутрішньої і зовнішньої динаміки. Дослідження динамічних 
процесів, з метою спрощення, в основному пов'язане з розглядом динаміки односту- 
пінчастих передач. Як у випадку одноступінчастих, так 1 багатоступінчастих передач 
знання частот власних коливань є визначальним для динамічного аналізу механічної 
системи. 

Механічні коливальні системи у класичній постановці складаються з типових ла- 
нок, які можуть бути представлені пружинами, демпферами 1 масоінерційними елемен- 
тами. Для розкриття фізичної сутності перебігу вібраційного процесу 1 отримання спект- 
рів вібрації використовують різні підходи, що базуються на розробці і дослідженні 
математичних моделей дискретних і континуальних систем |1|, |2). В залежності від 
характеру зовнішнього навантаження 1 пов'язаних з ним збурень, динамічні моделі по- 
діляються на детерміновані, стохастичні й евристичні (1, 21, І3). У детермінованих 
моделях зовнішні впливи 1 параметри розглядаються як невипадкові величини 1 функ- 
ції. У стохастичних моделях ці характеристики розглядаються як випадкові процеси, 
функції випадкових аргументів і випадкові величини. Керовані пружні властивості 
опор, в залежності від зовнішніх збурень, відображаються у евристичних моделях. 

За адекватністю опису динамічними моделями реальних процесів, на відміну від 
дискретних моделей, у яких маси, жорсткості та інші параметри розглядаються як зосе- 
реджені, в континуальних системах використовуються розподілені параметри. У змішаних 
моделях розподілені і зосереджені параметри поєднуються одночасно. Рух континуаль- 
них систем, як систем з розподіленими параметрами, описується диференціальними рів- 
няннями у частинних похідних, розв'язування яких для найбільш повних динамічних 
моделей зубчастих передач є складним, навіть у випадку застосування сучасних мето- 
дів комп'ютерного моделювання. У ряді випадків для визначення частот коливань і 
динамічних навантажень континуальних систем використовують алгоритмічні мето- 
ди кінцевих різниць і кінцевих елементів |3|. При цьому масо-жорсткісні, демпфуючі 
та інші параметри системи вважаються постійними, а їх зв'язки здійснюються у вузлах 
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області розбиття. Приймаючи переміщення вузлів за узагальнені координати системи, 
переміщення точок всередині елемента виражаються через переміщення його вузлів. 
Зазвичай, для континуальних моделей при моделюванні на мікрорівні, система розв'язу- 
ваних рівнянь руху має високий порядок, однак наявність великої кількості нульових 
елементів матриць жорсткості та інерції дещо спрощує розв'язок. Врахування пружних 
властивостей, фактора затухання коливань, наявності зазорів, нелінійності характерис- 
тик значно ускладнює динамічні моделі. 

При складенні динамічних моделей формується розрахункова схема, визначаю- 
ться місця і величини прикладених зовнішніх впливів (силових чи кінематичних), скла- 
даються рівняння руху динамічної системи у вигляді математичної моделі і обирається 
метод розв'язку або здійснюється підготовка початкових даних для комп'ютерного 
моделювання |3. 

Загальна методика розв'язування задач динаміки об'єктів зі стаціонарними зв'яз- 
ками грунтується на теоремах аналітичної механіки. Фундаментальним положенням 
є поняття принципу можливих переміщень, або принципу Лагранжа |4|. Об'єднуючи 
цей принцип з принципом Даламбера, Лагранж отримав загальні рівняння динаміки, з 
яких випливають основні диференціальні рівняння руху матеріальної системи 1 основні 
теореми динаміки. У (4) показано, що якщо рух голономної системи описується уза- 
гальненими координатами 4,45... 9, 1 узагальненими швидкостями 4; 45 Й,» ТО рів- 


няння руху мають вигляд: 


В райо прідасн (1) 
зчучва; 1ода, 


де Т - кінетична енергія системи, ПІ - потенціальна енергія, 4; - узагальнені координати, 
4; - узагальнені швидкості. Ці рівняння називаються рівняннями Лагранжа другого 


роду. Якщо рух відбувається у потенціальному полі і вважається, що П не залежить 
від узагальнених швидкостей 4, , то рівняння Лагранжа (1) набудуть вигляду |4Ї. 


4 д г аб бі 
топ перен а іаній» а 2 
и сн аа, ду й 0) 


/ 


і 

де ІА Т-П - кінетичний потенціал, або функція Лагранжа, п - кількість узагальнених 
координат. Якщо крім активних сил, що визначаються функцією П, на елементи системи 
діють сили, які цією функцією не можуть бути визначені, наприклад, сили опору різ- 
ного фізичного походження, то рівняння (1) можна подати у вигляді 


п 0, /-12,п. (3) 


У правій частині цих рівнянь знаходяться сили, що не визначаються функцією П. 

Перевагою рівнянь Лагранжа порівняно з іншими методами є те, що вони дають 
змогу розробляти математичні моделі динаміки за єдиною методикою. Переміщення 
елементів системи можуть бути як лінійними, так і кутовими, тобто узагальненими ко- 
ординатами системи можуть виступати як лінійні переміщення, так 1 кути повороту. 
Відповідні узагальнені сили - це звичайні сили і крутні моменти. Безпосередній аналіз 
зв'язку переміщень і сил на основі законів динаміки (Ньютона) і виведення відповід- 
них рівнянь руху пов'язаний з певними труднощами через векторний характер пере- 
міщень і зусиль. Тому для виведення рівнянь руху вважається доцільним використо- 
вувати рівняння Лагранжа другого роду у потенціальному полі |4|. Оскільки величини 
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Т 1 П є скалярними, то виведення рівнянь руху є досить простим і зводиться до зви- 
чайних операцій диференціювання. Необхідно відмітити, однак, що у простих випадках 
для виведення рівнянь руху можна безпосередньо використовувати закони Ньютона. 

Аналіз наукових досліджень. Актуальність дослідження динаміки зубчастих 
передач підтверджується роботами багатьох вітчизняних та зарубіжних дослідників. 
В основоположних роботах (5, с. 175, |6| були виявлені фази кромочного 1 серединно- 
го контакту зубців, вплив циклічної похибки, виконана оцінка динамічних навантажень 
у зачепленні, отримані спектри коливань зубчастої передачі. У роботах інших авторів 
динамічна модель була доповнена з урахуванням впливу пружної піддатливості опор; 
розглянуто частотний спектр параметричних коливань в залежності від змінної жорст- 
кості зачеплення 1 впливу масляної плівки; у окремих роботах відмічається суттєвий 
вплив на динамічну навантаженість і частотний спектр коливань пружних опор зубчастих 
передач. 

Метою даної статті є розробка просторової математичної моделі власних частот 
і форм коливань механічної системи класу циліндричних одноступінчастих евольвент- 
них зубчастих передач, з метою її подальшої комп'ютерної реалізації. Просторовість 
моделі забезпечується врахуванням шести ступенів вільності для зубчастих коліс 1 
трьох ступенів для кожного валу. 

Основна частина. Розробка математичної моделі здійснюється на основі попе- 
редньо створеної моделі динаміки, спрощена кінематична схема якої зображена на рис. 1. 


Рисунок 1 - Спрощена кінематична схема зубчастої передачі 


На схемі позначені: 1, 6 - приєднані маси двигуна і робочого механізму, 2, 5 - ді- 
лянки валів між приєднаними масами і зубчастими колесами, 3, 4 - зубчасті колеса. 

Запропонована попередньо модель динаміки дає змогу дослідити вплив мас зуб- 
частих коліс, пружність 1 довжину валів, на значення частот і форм власних коливань 
системи з урахуванням опору середовища, обумовленого наявністю тертя і демпфуючих 
властивостей масляної плівки у зачепленнях. За узагальнені координати динамічної 
моделі взято: 

Ф, Фі, Ф», фз - кути повороту приєднаних мас та ділянок валів навколо осей х/, х2; 

ф/.Ф/ Ф.Ф» - кути повороту шестерні та колеса навколо осей у),, у2 
Та 21,22 

хі У 2 Х», Уз 22 - лінійні переміщення шестерні та колеса вздовж осей х, у, 2. 

Конструктивні параметри механічної системи: 

7, Д- моменти інерції приєднаних мас двигуна І та навантаження 2 відносно 
осі х (рис. 1); 

з,» - моменти інерції відносно осі х ділянок валів 2, 5 між зубчастими коле- 
сами 3, 4 та приєднаними масами |, 6 відповідно (рис. 1); 


яд» г 7 мОМенти інерції зубчастих коліс 3, 4 відносно осі у (рис. 1); 


. г - МОМенти інерції зубчастих коліс 3, 4 відносно осі 2 (рис. 1); 

ті, т - маси зубчастих коліс 3, 4 (рис. 1); 

М, Мз - крутний момент двигуна І та момент опору механізму навантаження 6 
відповідно; 
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Ко, Ко, П, Із - крутильні жорсткості та довжини валів 2, 5 (рис. 1); 
СІ, С, - жорсткості опор вхідного вала 2 по осі 2 (рис. 1); 


СІ, С, - жорсткості опор вхідного вала 2 по осі у (рис. 1); 


С. - жорсткість опори вхідного вала 2 по осі х (рис. 1); 


27» Су, - жорсткості опор вихідного вала 5 по осі 2 (рис. 1); 


' 


зу» С, - ЖОрсткості опор вхідного вала 5 по осі у (рис. 1); 


С,» - жорсткість опори вхідного вала 5 по осі х (рис. 1). 


Таким чином, всі конструктивні, геометричні та динамічні параметри відобра- 
жені практично повною мірою. Відповідно, вектор узагальнених координат механічної 
системи містить сукупність кутових та лінійних переміщень і має вигляд: 4 - (Ф, ФІ, Ф», 


Фу, фіз фоз Хр Різ Зо УіВЯ : 

Отже, вектор узагальнених координат містить чотирнадцять компонент, які опи- 
сують крутильні і лінійні (поздовжньо-поперечні та згинальні) коливання елементів 
механічної системи. 


Методика складання рівнянь Лагранжа (1) - (3) містить такі кроки: 
-2 


4 


1. Визначення кінетичної енергії системи Т « У , де Т. «кінетична енер- 


І 
гія окремих елементів системи. 
2. Визначення потенціальної енергії системи П- РІ , де П.,- потенціальна 


І 
енергія окремих елементів системи. 
3. Знаходження частинних похідних від кінетичної енергії по узагальнених швид- 


Об. пори туни а|оТ божа 
костях з і повної похідної за часом РОЙ Кількість узагальнених координат 
Ч; Ї 4; 


відповідає кількості ступенів вільності системи. 
4. Визначення частинних похідних від кінетичної енергії по узагальнених коор- 


динатах | --- 
94; 
5. Визначення узагальненої сили О), яка для консервативної системи визначиться 
2 2 | 3 б ща Е ОП 
як частинна, похідна від потенціальної енергії по узагальненій координаті 0, « -| - -- 
04; 
Виконаємо процес створення математичної моделі по кроках. 
1. Кінетична енергія системи запишеться у вигляді: 
ЕР зго РАВ УЛ зум2 292 42 
2 фут (Ху ВУ ка, ун уз/щфУчУ/ьфФ Я 
(4) 


222142 зум2 і292 -2 42 

ет, (У 25) фо Фо ФЛ фу. 

За відсутності приєднаних мас, тобто у випадку ізольованої зубчастої передачі, 
члени ./Ф/;./Ф; будуть відсутні. За відсутності осьових переміщень, тобто без ураху- 


: : є .2 .2 
вання поздовжніх осьових сил, відсутніми будуть члени т,Х; та т,Х.. 


2. Визначаємо потенціальну енергію системи: 
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2пе-К(Ф фу) Кк (Ф» фу) СІ, (2 Іф - СІ, (2 чІф/) за ОГО -1фІ) ї 


ро СІ, (2, 4 1фІ) ж (У ІФ) СІ, (и я фу С; (з -1фу) нСІ, (у» фі) на 
1 г укр , 
нео о Каїчядунів офі ять зтрарь НЕВМОх зандрнічи 0 9) 
1 Б З 2 
ни з узОЇ фі тив | жд) Р С, 


віп о, 


3. Проведемо операції диференціювання згідно з формулою (1). Похідні від кіне- 
тичної енергії: 


б з пічоо ПРА А п - по обі 
зо у- гав 9 - за у- фіз? -7афі 
х б у-ту У 2. уз тубі У 9 уз тур Я 551 -т,), 


4. Похідні від кінетичної енергії по узагальнених координатах у всіх випадках до- 


ОТ 
рівнюють нулю, тобто -- - 0. 
4; 
5. Проведемо операції диференціювання від потенціальної енергії. Відповідні час- 
тинні похідні набудуть вигляду: 


ОП дП ДОП 
Рон 102 Ф3 чі к (Ф» Фа); 


ОП 
аЗК(Ф тео В 2 (Ф, - фі) Я С, (Гу, ХВ; 


дФ 93 дФ, 2 
ОП і й ії ; 
29 -ос,2,з Пон с,єі з іФіЛ,чсС,() ха х Ва; 
і 
ОП ' " у 
до - ос, 19ІЛ, с, І9ФіЛ, - С() х с х Ву; 
2 
оп '" 2 "її І 
до жа або сейфу Да З су: зд ОО о Ро ке Р о 
1 
ОП ' 52 б Ж 
29: неба СаО ЛОРОТЛЕ с, (У» Т І.Ф, )1, - С,(п)х пу, х 7 В, х Кб; 
2 
о зс,х я С (Фає в, -1)х Ро СЯ ее хо-СбФІЄВ.ФІ1)х В. 
Ох, дх; 
ОП ' щі "ї І 
ду з сії 7 ФІН ск 1ФІ)чН С, х В; 
1 
ОП ' Р " І 
су - с5(У» -ЛФі)ч су» - І,ф;)ч СО) х Ку; 
2 
ОП ' у и 
да - са; Фі сце чіф/)-сС(0хк,; 
1 
ОП ' "І у 
з - с,(2, - 93) су (2, 1 1Ф2)- С (0) х В,,; 
2 
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Коефіцієнти КБ; - К;г являють собою циклічно повторювані функції. На підставі 
формули (1) 1 проведених операцій диференціювання, прирівнявши ліві і праві частини, 
отримуємо систему 14 диференціальних рівнянь руху динамічної системи зубчастої 
передачі, без врахування фактора демпфування системи. Отримана математична мо- 
дель має вигляд: 
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Кожне з рівнянь системи окремо описує: 

- рівняння І (відносно ф ) - крутильні коливання ділянки первинного вала 2 між 
приєднаною масою двигуна | та зубчастим колесом 3 (рис. 1); 

- рівняння 2 (відносно Ф, ) - крутильні коливання колеса 3 навколо осі х; (рис. 1); 

- рівняння 3 (відносно 9, ) - крутильні коливання колеса 4 навколо осі х» (рис. 1); 

- рівняння 4 (відносно о, ) - крутильні коливання ділянки вторинного вала 5 між 
приєднаною масою навантаження 6 та зубчастим колесом 4 (рис. 1); 

- рівняння 5 (відносно ф/ ) - крутильні коливання колеса 3 навколо осі у; (рис. 1); 

- рівняння 6 (відносно ф» ) - крутильні коливання колеса 4 навколо осі у» (рис. 1); 

- рівняння 7 (відносно ф )- крутильні коливання колеса 3 навколо осі 7, (рис. 1); 


- рівняння 8 (відносно ф, ) - крутильні коливання колеса 4 навколо осі 2. (рис. 1); 
- рівняння 9 (відносно х/) - поздовжні коливання колеса 3 вздовж осі х; (рис. 1); 
- рівняння 10 (відносно х2) - поздовжні коливання колеса 4 вздовж осі х» (рис. 1); 
- рівняння 11 (відносно у,) - поперечні коливання колеса 3 вздовж осі у; (рис. 1); 
- рівняння 12 (відносно у) - поперечні коливання колеса 4 вздовж осі у» (рис. 1); 
- рівняння 13 (відносно 2/) - поперечні коливання колеса 3 вздовж осі 7; (рис. 1); 
- рівняння 14 (відносно 22) - поперечні коливання колеса 4 вздовж осі 2» (рис. 1); 
Таким чином, перші 8 рівнянь системи описують крутильні коливання елемен- 
тів механічної системи, решта 6 - поздовжньо-поперечні коливання зубчастих коліс. 
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Висновки 


На основі попередньо розробленої моделі динаміки, з використанням рівняння 
Лагранжа другого роду, створена просторова математична модель (6) власних частот 
та форм коливань механічної системи класу циліндричних, одноступінчастих, евольвент- 
них, косозубих зубчастих передач. Запропонована математична модель являє собою 
систему 14 диференціальних рівнянь другого порядку, відносно узагальнених координат 
динамічної моделі. В залежності від застосованого алгоритму розв'язок даної системи 
диференціальних рівнянь може бути отриманий у вигляді амплітудно-частотного спектра 
частот та форм власних коливань динамічної системи (частотна область) та у вигляді 
зміни у часі кожної з узагальнених координат (часова область). 

Для першого випадку задача пошуку розв'язку системи зводиться до обчислення 
власних чисел характеристичної матриці системи рівнянь (спектра матриці). З цією ме- 
тою можуть бути застосовані відомі чисельні методи, наприклад, метод А.М. Данилев- 
ського, Леверрьє - Фаддєєва, невизначених коефіцієнтів тощо. 

У другому випадку для отримання розв'язку необхідно скористатись методами 
покрокового інтегрування типу Рунге, Кутта - Мерсона, Ейлера, тобто звести рішення 
рівнянь системи (6) до рішення задачі Коші. Це передбачає заміну системи з 14 дифе- 
ренціальних рівнянь 2-го порядку системою з 28 рівнянь 1-го порядку і визначення 
початкових умов, зазвичай з умови статичної рівноваги описуваної механічної системи. 
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П.В. Дяченко 

Пространственная математическая модель собственньгх частот и форм колебаний 
механической системьі класса одноступенчатьх, звольвентньтх зубчатьтх передач 

В статье предлагаєется пространственная математическая модель динамики механической системь класса 
одноступенчатьх, цилиндрических, косозубьтх, звольвентньх зубчатьїх передач. Предлагаємая модель 
представляєт собой систему четьшрнадцати дифференциальньх уравнений второго порядка относительно 
обобщенньгх координат, составленную на оснований уравнений Лагранжа второго рода, которь описьтвают 
пространственнье колебания зубчатьїх колес механической системьі. 


Р.У. Буаспепко 
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